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^ ■ Abstract 

>>: 

The QCD factorization (QCDF) method usually contains infrared divergences which introduce 



large model dependence to its predictions on charmless B decays. The amplitudes of charmless B 
decays can be decomposed into "tree" and "penguin" parts which are conventionally defined, not 
from the topology of the dominant diagrams, but through their associated CKM factors V* b V uq 
Ph| and V t * b Vt q , respectively, with q = d,s. We find that for B^ s — > ir^K^ decays, the "tree" amplitude 



can be well estimated in QCDF with small errors, as the endpoint singularities have been canceled 



to a large extent. With this as the only input from QCDF and combined with flavor symmetry, the 
branching ratio of B s — > ir + K~ is estimated to be significantly larger than the CDF measurement. 



This contradiction could be solved if the form factor F BsK is smaller than the light cone sum rules 



in 

CN ' estimation or the "tree" amplitude has been over estimated in QCDF. The latter possibility could 

o ; 

happen if charming penguins are nonperturbative and not small, as argued in soft collinear effective 
theory. To differentiate between these two possibilities, we examine the similar B s — > p + K~ decay 
with the same technique. It is found that a large part of the uncertainties are canceled in the ratio 
B(B S -»■ p + K-)/B(B s -> 7T+K-). In QCDF, it is predicted to be 2.5±0.2 which is independent on 
the form factor. However if charming penguins are important, this ratio could be very different from 
the QCDF prediction. Therefore the ratio of these two branching ratios could be an interesting 
indicator of the role of charming penguins in charmless B decays. 
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I. INTRODUCTION 



Charmless hadronic B decays provide various possibilities for the determination of the 
CKM angles, but often the efforts are hindered by the low energy QCD dynamics which may 
lead to large, sometimes even uncontrolled, systematic uncertainties. To proceed further, 
one may resort to either factorization methods or flavor symmetries. However it is known 
that here factorization only holds to the lowest order of power expansion in l/m^, which 
means there is no (known) model-independent way to estimate power corrections in l/mj, 
systematically. For flavor symmetries, in general SU(3) flavor symmetry breaking effects may 
also introduce large uncertainties though isospin is a good symmetry at the tew percent level. 
Recently, Descotes-Genon et al. [l( observed that, by the combined use of QCD factorization 
(QCDF) and flavor symmetries, it is possible to obtain new sum rules for Bd yS — > KK decays 
with better control on systematic uncertainties. This, together with the first observation 
of B s — > tt + K~ by CDF Collaboration [2], stimulates us to investigate Bj — > n^K^ + and 
B s — > n + (p + )K~ decays in the spirit of Ref. 

Factorization aims to separate the calculable short distance effects from the nonpertur- 
bative long distance part. In addition, the nonperturbative structure can be simplified as 
form factors, decay constants and light-cone distribution amplitudes (LCDAs) in hadronic 
B decays at leading power expansion in I/mi,. For charmless B decays, there are three fac- 
torization approaches available in the market: QCDF |^-5], the perturbative QCD method 
(PQCD) 6-8] and soft collinear effective theory (SCET) [9-11]. We shall not discuss PQCD 
in this paper. For QCDF, phenomenologically power suppressed corrections such as chi- 
rally enhanced power corrections and annihilation topology have to be included to fit the 
experimental data. Unfortunately these terms contain endpoint singularities which break 
factorization and can only be estimated in a model-dependent way. In SCET, the chirally 
enhanced power corrections and weak annihilation diagrams are calculable jl^ . [ijj], while 
the diagrams with internal charm quark loops are claimed to be nonperturbative. Con- 
sequently these charming penguin amplitudes 1 can only be determined by fitting to the 
data. Generally as weak annihilation and charming penguins have the same topology, it is 
hard in practice to tell which one is really important. In this paper we find that the ratio 



1 It was first proposed in 15l-ll 7T| that the charming penguins may play an important role in charmless B 
decays. 
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B(B S — > p + K~)/B(B s — > 7i + K~) may provide some insight into this issue. 

This paper is organized as follows. In the next section, we first discuss B d)S — » tt t K ± 
decays with a minimal use of QCDF combined with flavor symmetry. We then extend 
the discussions to B d — > n~K* + and B s — > p + K~ decays in the third section. Finally we 
conclude with a summary in section IV. 



II. B d>s -> ir^^ DECAYS 

Recently the CDF Collaboration reported new results on the branching ratios and direct 



CP violation for B® s decay channels {2], among them is the first observation of B s — > n + K 
decay: 

B(B S 7v + K-) = (5.0 ± 0.7 ± 0.8) x 10~ 6 , (1) 
Acp(B s — > n^K~) = 0.39 i 0.15 i 0.08 , (2) 



where the first errors are statistical and the second systematic. It was first proposed in 12] 
to determine the CKM angle 7 from B^ s —> tt t K ± decays using the U-spin symmetry. It was 
further noticed in that here the flavor symmetry breaking effects should be unusually 
small since the strong phases from final state interactions are exactly the same due to the 
charge conjugation symmetry of the final states. Therefore it could be a robust test of the 
Standard Model vs New Physics to check the equality of the direct CP asymmetries of these 
two decay channels. In this paper we will instead focus on their branching ratios, together 
with (Bd — > tt~K* + , B s — > p + K~) decays, with a minimal use of QCDF combined with 
flavor symmetry 2 . 

For tiK channels, the decay amplitudes can usually be expressed as 

A(B d n-K + ) = At K (V: b V us T d + V; b V cs P d ) , (3) 
A(B S -> n + K-) = Al K {V: b V ud T s + V; b V cd P s ) , (4) 



with 



2 A recent analysis of B c i. s —> Kir decays using flavor symmetry can be found in 19 1 
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T d,s and P d ' s j multiplied by a common factor A^, denote the "tree" and "penguin" am- 
plitudes accompanied by the CKM factors V* b V uq and V* b V cq , respectively, with q = d,s. In 
QCDF fl, 

T = ai + at + a\ ew + ft - \^ ew , 

P = at + al ew + (3 C 3 - l -Pl ew , (5) 

where the weak annihilation contributions are contained in j3 terms, while all the other parts 
are grouped into „ terms, ineluding the vertex eorreetions and the hard speotator mattering 
contributions. The explicit expressions of a s and p a can be found in [5J. In QCDF, 
the endpoint singularities appear in annihilation topology and hard spectator scattering 
diagrams, which introduce significant uncertainties to a's and /3's. Therefore neither T nor 
P can be predicted reliably in QCDF with small errors. However, by implementing the 
unitarity condition of the CKM matrix, the above decay amplitudes can be reexpressed as 

A(B d -> n-K + ) = At K (V: b V us f d - V: b V ts P d ) , (6) 
A(B S n + K-) = A s wK (V* b V ud f s - V*V td P s ) , (7) 

with T = T — P. It was first observed in |l| that for B d s — > K°K° decays, although both 
T and P contain endpoint singularities, T is infrared safe as the endpoint singularities in T 
and P cancel completely. Therefore T can be estimated reliably in QCDF. However such 
case is rare in B decays. For B — > K7i decays, using Eq. (jSJ), one obtains 

Interestingly, for every bracket in the above equation, the endpoint singularities are canceled 
exactly, so the only residual endpoint singularity lies in ol\ which comes from the hard 
spectator scattering diagrams. Therefore, although T in B — > Kit decays is not totally 
free of infrared divergence, in practice it could be estimated in QCDF with much smaller 
uncertainty. 

Numerically, the errors of T are dominated by the residual endpoint divergence from hard 
spectator scattering diagrams which is usually modeled as 



(9) 
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and the parameter Xb, which parameterizes the integral of the B meson wave function as 

f 1% ^ ? • <"» 

JO i — ? 

To estimate the uncertainties, we will take Xb = 300 ± 100 MeV, < pn < 1, < < 27r 
and adopt the default values of [^j for all the other parameters. We then find 

\f s \ = 0.99i°;°5 , \T d \ = l-00i°;°4 , (11) 

which respects the U-spin symmetry T s = T d at one percent level. Thus in the following we 
will not differentiate between T s and T d and simply take in the flavor symmetry limit 



-id\ 



\T S \ = T = 0.99±g;gi . (12) 



It is worth reminding that charming penguins, accompanied by the CKM factors V* b V cq , will 
contribute to both T and P. So if charming penguins are nonperturbative and not small as 
argued in SCET, the above estimation of T could change. 

As mentioned earlier, in QCDF the "penguin" amplitude P is significantly affected by 
the endpoint divergences and therefore hard to be estimated in a reliable way. So the above 
equation of T is the only result that we will adopt from QCDF. In this sense B^^ — > Kir 
decays are analyzed with "a minimal use of QCDF". For all the rest we will resort to flavor 
symmetry. As observed in [3], the flavor symmetry breaking effects in B^s 7r T K ± 
should be unusually small due to charge conjugation which is a strict symmetry for the 
strong phases from final state interactions. So by taking the factors A^ S K out which have 
included explicitly part of the U-spin symmetry breaking effects, it should be reasonable 
to take the approximation of the U-spin symmetry limit P d ~ P s = P. This amounts to 
neglect hard spectator scattering contributions sensitive to the difference between B^ and 
B s distributions, and also annihilation contributions when the gluon is emitted from the 
spectator quark. Then the decay amplitudes can be rewritten as 

A(B d -> n-K + ) = At K (V u e^\f\ - V t \P\e lS ) , (13) 
A(B S ->■ W+K-) = K K (yy^\f\ + Vie^\P\e a ) , (14) 

where V q = |V^* b V^ s | and V' q = \V* b V q d\ for q — u, t. 5 is the relative strong phase between P 
and T: 

P/f = -\P/f\e iS , 



where the convention is chosen such that 5 = in the naive factorization limit. It is then 
straightforward to get a well-known relation between the direct CP violations of Bd. s — > kK 
decays 

A CP (B S -> K+K-) _ B(B° -> n-K+) t(B s )(A^ k ) 2 
A cp (B°^it-K+) B(B s ^n+K-) T (B°){Ai K y 1 ' 

Concerning the currently large experimental uncertainty of Acp{B s — y n + K~), this relation 

is well consistent within errors with the current data of B s — y ttK in Eqs.(l,2) and of 



B r] -> ttK 



20] 



B(B d ->■ 7v-R + ) = (19.4 ± 0.6) x 10~ b , A CP (B d -> 7r"K + ) = (-9.81^)% . (16) 

In this sense, there are three independent experimental observables: two branching ratios 
and one direct CP violation. If the related form factors are taken as input, there are also 
three unknowns: CKM angle 7, the "penguin" amplitude |P|, and the strong phase S. Then 
in principle, it should be possible to extract 7 from B d s — > ir T K ± decays, if the QCDF 
knowledge of T is reliable. However in practice, it is hard to obtain a strong constraint 
on the CKM angle 7 through this way, especially considering the uncertainties of the input 
parameters and also the mild dependence of the branching ratios on the angle 7. As currently 
no clear evidence of new physics beyond SM exists in rare B decays, we will instead take the 



global fit determination of 7 = (67.9^38)° from CKMfitter 2l| as input. Then our focus is 
what information one may obtain about the soft QCD dynamics in B — y Kir decays. 

Now only the "penguin" amplitude |P| and the strong phase 5 are unknown, which can 
be extracted solely from the observables of B& — > tt~K + decay in Eq. ( JT6l) . With the CKM 
parameters taken from CKMfitter [2l| : 

A = 0.2251 , \V ts \ = 0.040421™ , |K 6 | = 0.00351+™ , (17) 

, the determination of \V ub \F BlT = (9.1 ± 0.7) x 10~ 4 |22j from the semi-leptonic B — y ir£u 



decay and r Bd = 1.525 ± 0.009 ps |20|, we obtain 



|P| = 0.129 ±0.012 , S = -(18.9 ±2.9)° , (18) 

with the uncertainties coming mainly from the variation of |\4b|P B7r . Notice that the deter- 
mination of |P| is rather insensitive to the change of T. For instance, about 20% change of 
|T| from 0.99 to 0.80, with all the other parameters unchanged, will lead to only 1% varia- 
tion of the central value of |P| from 0.129 to 0.128. This is because, as a penguin dominant 



io 2 |^| 


u 


ltf\V ub \F B « 




lO 6 ^^ -»• TT+if ) 


4 n42+ 037 


131 MeV 


9.1 ±0.7 


1 472 +0.024 
i - 4 ' z -0.026 P s 


23.1 ± 1.0 


io 3 |^| 


fK 


pB a K jpBir 






8.59«;1 


160 MeV 


i ' io -0.09 


1.638 ±0.011 ps 





TABLE I: Input parameters in numerical calculations of B s — > ir + K decay. The CKM parameters 
are taken from CKMfitter 21[ , The 5-meson lifetimes and the branching ratio are taken from [2j 



and the form factors are given in 



22, 



23]. 



decay, the branching ratio of Bd — > tt~K + decay largely determines the magnitude of P, 
irrespective of how much the charming penguins may contribute to P. 

It is also known that B + — > tt + K° decay is almost pure penguin process, in addition its 
penguin amplitude is equal to that of B d — > n~K + decay except for the negligible color- 
suppressed electroweak penguin contributions. As a consistency check, we do find that the 
penguin amplitude of B + — > tt + K° is \P\ = 0.130 ± 0.012, which is in well agreement with 
Eq. (11 8p . Therefore the above estimation of the "penguin" amplitude \P\ in Eq. ( !T8|) should 
be, up to few percent, independent on whether the charming penguins are important or not. 

Taking Eq. (j!8p as input, we can now estimate the branching ratio of B s — > ii + K~ decay. 
The CP-averaged amplitude square can be expressed as 



\A(B S 



tt+K' 



K 



) 2 W'T\ 2 



1 



Vl_P_ 



+ 2 



V[P 



cos 5 cos(/3 + 7) 



(19) 



Due to the accidental fact that + 7 ~ 90°, the "tree" -"penguin" interference term in 
B s — > ir + K~ decay is negligibly small. Thus approximately the strong phase 5 does not 
affect the branching ratio and one has 



\A{B S vr + iv-)| 2 ~ (Al K f ( |KT| 2 + |F/P' 2 



(20) 



If one observes that the \P\ 2 term can be directly related to the pure penguin process 
B+ n + K° as 

(21) 



numerically one can then obtain 



B{B S vr+iT 



F B a K ^ 2 
pBn 



(6.2 ± 1.0)|T| 2 + 0.6 



x 10" 



(22) 
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with the input parameter listed in Table I. Here the uncertainty of \P\ term, which can be 
estimated by using Eq. (|2T|) . is much smaller than that of T term and has been therefore 
neglected. We want to stress that Eq. (I2"2"j) is not just valid in QCDF, but rather be a model 
independent result to a large extent. 

The form factors ratio has been estimated in light-cone sum rules [23] , as shown in Table 
I. As discussed earlier, T can also be estimated in a relatively reliable way in QCDF in Eq. 
f JTTj) . then we get 3 

B th (B s -> 7t + K-) = 8.8±i;t±?:i x 10" 6 , (23) 

where the first error comes mainly from the uncertainties of jV^I-F^ 71 " and T, while the 
second error comes from the form factors ratio. This estimation is significantly larger than 
the CDF measurement (5.0 ± 1.1) x 1CT 6 . 

Notice that |Ki&|-F 1 ' B7r is extracted from B — > tt£v decays in a model independent way. 
Thus if the updated experimental measurement does not change greatly, it must be that 
either the form factors ratio is much smaller than the light-cone sum rules prediction or T is 



overestimated in QCDF, or both. The possibility of a smaller ratio of the form factors 



nas 



19| and also in a recent QCDF analysis 



24 



25| 



been discussed in a flavor symmetry analysis 
where the central values of F BaK = 0.24 and F Bn = 0.25 are taken. The latter possibility of a 
smaller T can not be easily realized in QCDF, though it could happen in SCET if charming 
penguins are not small. In practice, however, it is difficult to tell experimentally which 
possibility is chosen by nature. This is because T is not directly observable in charmless B 
decays. For the form factor F BsK , in principle it could be determined via semileptonic B s — > 
Kiv decays. But this measurement is not easy at hadron colliders, such as the Tevatron and 
the LHC. For the future super B factories, the challenge is to accumulate enough B s mesons 
for precise measurement of the kaon momentum spectrum of this semileptonic decay. F BsK 
may also be predicted by lattice QCD computations at small recoil, while unfortunately 
charmless B decays happen at large recoil. 



3 A direct QCDF calculation in |5j gives 10.2tf;f±f;|±?;£t2'.7 x 1CT 6 as their "default results". Here the 
first error corresponds to the variations of the CKM factors and the second error comes mainly from the 
variation of the form factor. The readers are referred to [5( for more details of the analysis. 
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III. B d -> tt-K*+ AND B s -> p+if- DECAYS 



To find a way out of this problem, we may examine similar decay channels Bd — > n~K* + , 
B s — > p + K~ and B + — > tt + K*°. At first sight, it may seem better to discuss B^^ — > p T K ± 
decays as the final states are related by charge conjugation symmetry. However these two 
decay channels are actually not related by U-spin symmetry. Again we may express the 
relevant decay amplitudes in the following form: 

A(B d -> n-K* + ) = A« K *{Vj"i\f K *\ - V t \P K *\e iSK *) , (24) 
A(B S -> p+K-) = A pK {Vy\f p \ + V;e~^\P p \e i5 ") , (25) 



with 



A nK * = V2G F m BdPc f K «F B7T , A p k = V2G F m BsPc f p F 



B S K 



Here p c is the cm. momentum of the final state mesons. In QCDF, the expression of T p is 
formally the same as Eq. ([8]) of T, though the explicit formulae of the coefficients a's and 
/3's are different which can also be found in [5j. Since now the final states irK* and pK are 
not related to each other under charge conjugation, there is no reason to expect that here 
the U-spin symmetry breaking effects would be particularly small. This means the difference 
between (T p , P p , 5 P ) and (Tk*, Pr*, 8k*) could reach the level of 20% and can not be simply 
ignored. But if we focus only on the B s — > p + K~ decay, the strong phase S p again does not 
affect the CP-averaged branching ratio 

\A(B S p+K-)\* ~ (A s pK ) 2 (|KT/ + \V;P P \ 2 ) (26) 

in good approximation due to the accidental fact that (3 + 7 ~ 90°, just like the case of 
B s — > n + K~ decays. Similarly we find \T p \ = 0.99^qq5 with the uncertainties mainly from 
Xb and the hard spectator parameters p^ v and 4>h V ■ We will soon find the ratio T p /T to be 
useful. As they share the common parameter A#, their errors are partly correlated and the 
uncertainty related to Xb will be canceled in the ratio. But the hard spectator parameters 
are independent for B s — > PP and B s — > PV decays. Concerning this, we find the ratio to 
be 

T 

p - 1.00 ±0.04 



T 

with the uncertainty coming dominantly from the variation of hard spectator parameters. 
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Then the last missing piece required is the "penguin" amplitude \P P \. Notice that B s — > 
p + K~ decays are dominated by the "tree" amplitude T p . Furthermore it is known that 
generally the penguin amplitudes of B — > PV decays are smaller than those of B — > PP 
decays. As the \P\ term contribute about 10% of the final branching ratio of B s — > -k + K~ 
decay, it is expected that the contribution of \P p \ term should be even smaller. Therefore for 
such a small term, an estimation based on SU(3) flavor symmetry should be enough which 
will not introduce large uncertainty to the final branching ratio result. 

In the flavor symmetry limit, one has \P P \ = \Pk*\- Then one may estimate \Pk*\ from 
either B d — > n~ K* + decay or the almost pure penguin B + — > tt + K*° decay, just like the 
case of B — > ttK decays. Taking f p = 216 MeV, fx* = 220 MeV and the parameters listed 



0.066i°;™6 using 



in Table I, we find \Pk* 

B(B d n-K* + ) = (8.6 ± 0.9) x 10 



20] 



Acp(B d -> ir~K* + ) = (-18 ± 7)% , (27) 



which is consistent with the determination of \Pr* 
(9.9±^;|) x 10" 6 [20]. Then we obtain 



n nfi4+ 007 



\A{B S -> p + K- 



(A s pK f [ KT P \ 2 + 



from B(B + n + K*° 



A 



ttK* 




(28) 



Numerically the corresponding branching ratio can then be estimated 



B{B S -> p + K~) 



pBn 



;i5.7±2.4)|T p | 2 + 0.4 



x 10" 



(29) 



where the uncertainty from \P p \ term is much smaller than that from T p term and has 
therefore been ignored. We want to stress again that Eq. (1291) . just like Eq. ( )22l) . is a model 
independent result to a large extent. 

Notice that B s — > 7i + (p + )K~ decays depend upon many common parameters such as 
F BaK , iKibjF' 677 etc., their uncertainties are highly correlated. It is then clear that the ratio 
of the branching ratios should have much smaller error. Actually the error of the ratio 
depends almost solely otlT p /T 



B(B S -> p+K- 
B(B S ->• 7T+K- 



f 2 

P 



T 



T 



1 + 0.4/(15.7|T p 



2.5 



T 



T 



(30) 



l + 0.6/(6.2|T| 2 ) 

In the above estimation, we have adopted flavor symmetries to relate the "penguin" ampli- 
tude of B s — > n + (p + )K~ to that of B d — > n~K^ + , which may introduce 20% level of uncer- 
tainty. But fortunately as tree-dominant decays, the ratio B(B S — > p + K~)/B(B s — > n + K~) 
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~ 2.5 ±0.1 ±0.2 , (31) 

QCDF 



is not sensitive to the variation of the relevant "penguin" amplitudes. As an illustration, we 
assign conservatively 30% uncertainty to the "penguin" amplitudes to take account of the 
flavor symmetry breaking effects, and adopt the QCDF estimation \T p /T\ = 1.00 ± 0.04, it 
is found that 

B(B S -» p+K- 
B(B S -»■ 7i+K- 

with the first error comes from the flavor symmetry breaking effects and the second error 
comes from the variation of \T p /T\. While in SCET, if the charming penguins are not small, 
it could alter T and T p . It is worth reminding that there is no known reason, at least to 
our knowledge, to expect that the charming penguins of B s — > n + K~ and p + K~ channels 
should be the same. This means the ratio B(B S — > p + K~)/B(B s — > tt + K~) in SCET could 
be quite different from E q . ( I3T1) . For example, a recent global fit analysis of charmless B 



decays using SCET gives 



26| 4 (in units of 10' 



B(B S -> TT+K-) = 5.7±S:i±8:| , B(B S p+K~) = 7.6l°;?i°J (Solution I) (32) 
B(B S 7T+K-) = 5.5t° 6 +° i , B(B S p+K~) = 10.2±g;f±g:| (Solution II) (33) 

In these two solutions, the ratio B(B S — > p + K~)/B(B s — > n + K~) equals to 1.3 or 1.9, 
significantly different from the QCDF prediction 2.5 ±0.2. Therefore a precise measurement 
of this ratio in the near future could be an interesting test of the importance of the charming 
penguins in charmless B decays. 



IV. SUMMARY 

In this paper we first studied B^s — > ti t K ± decays using QCDF combined with flavor 
symmetry. In general, the QCDF calculations contain infrared divergence when weak an- 
nihilation topology and hard spectator scattering diagrams are included. But if we express 
the decay amplitude as 

a(b s n + K-) = K K {v: b v ud f - V t * b V td P) , 

The "tree" amplitude T can be (relatively) reliably estimated in QCDF to be 0.99^q;q5. 
This is because the endpoint singularities inside T has been canceled to a large extent. In 

The results of tt + K~ were provided by Wei Wang, one author of [26(, in a private communication. 
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addition, the CP-averaged branching ratio is insensitive to the relative strong phase between 
T and P because the interference term is proportional to cos(/3 + 7) which is (accidently) 
very close to zero. Noticed that the final states of B^ s — > tt t K ± are invariant under charge 
conjugation, impl ying that the flavor symmetry breaking effect should be unusually small 
as pointed out in 18). So the penguin amplitude \P\ can be determined from — > ti~K + 
decay with small error. Then we obtain 



^B S K 



2 



B(B S ir + K~) = ( ) (6.2 ± 1.0)|T| 2 + 0.6 x l(T b . (34) 



With the light cone sum rules estimation [23j F BsK /F B7T = 1.15lg;^ and the QCDF pre- 
diction of T, the expected branching ratio would be 8.8^'g x 10 -6 , much larger than the 
CDF observation (5.0 ± 1.1) x 10~ 6 . One possibility is that the form factors ratio F Bs /F Bn 
has been overestimated by the light cone sum rules calculation. But it is also possible that 
the QCDF prediction of T is too large. This could happen, for example, if the charming 
penguins play an important role in charmless B decays. 

To differentiate between these two possibilities, we examined the similar decay channels 
B s — > p + K~, B d — > n~K* + and B + — > 7r + K*° with the same method. We find a large 
part of uncertainties has been canceled in the ratio B(B S — > p + K~)/B(B s — > n + K~) which 
equals approximately to be 2.5\T p /T\ 2 . In QCDF, \T p /T\ is also well predicted to be about 
one, so we get an interesting prediction B(B S — > p + K~)/B(B s — > n + K~) = 2.5 ± 0.1 ± 0.2 
in QCDF, with the first error reflects the flavor symmetry breaking effects and the second 
error reflects the uncertainty of \T P /T\. However if charming penguins are not small, this 
ratio could vary significantly from 2.5, which means the ratio of these two branching ratios 
could be useful to clarify the role of charming penguins in charmless B decays. 
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